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We develop a 7r-electron effective field theory (71--EFT) wherein the two-body Hamiltonian for 
a 7r-electron system is expressed in terms of three effective parameters: the 7r-orbital quadrupole 
moment, the on-site repulsion, and a dielectric constant. As a first application of this 7T-EFT, we 
develop a model of screening in molecular junctions based on image multipole moments, and use 
this to investigate the reduction of the HOMO-LUMO gap of benzene. Beyond this, we also use 
7T-EFT to calculate the 7r-electron contribution to the van der Waals interaction between benzene 
and a metallic electrode. 



Owing to the profound versatility of the carbon-carbon 
bond, organic molecules form the basis for a myriad of po- 
tential nanotechnology applications. Many of these make 
use of the ability of conjugated organic molecules to con- 
duct electricity, in which case the system of delocalized tt- 
electrons plays a role analogous to that of the conduction 
band in a conventional semiconductor. In such devices, 
the most important degrees of freedom from a technolog- 
ical perspective are those associated with these current- 
carrying 7r-electrons. Although the standard paradigm 
for transport calculations in molecular junctions-based 
on local approximations within density functional the- 
ory (DFT)-includes these degrees of freedom, it does 
so amongst many others. This comes at the expense of 
an exact treatment of intramolecular many-body interac- 
tions, which to date has made it challenging to capture 
important interaction effects such as Coulomb blockade 
within the context of DFT. 1 

An alternative approach is to formulate a model in- 
cluding only the degrees of freedom essential to describ- 
ing the 7r-electron dynamics, thereby reducing the over- 
head associated with an exact treatment of many-body 
interactions. This begins with the observation that pro- 
cesses in systems of 7r-electrons take place at characteris- 
tic length, energy and time scales all ultimately dictated 
by the strength of the 7r-electron bond. Intuitively, one 
expects that only degrees of freedom with scales com- 
parable to these need to be explicitly included. Semi- 
empirical models based on this notion have been in use 
for over fifty years j^r— and work to improve their accu- 
racy is ongoing However, since these are based on ad 
hoc parameterizations^r— of interparticle Coulomb in- 
teractions that do not satisfy Maxwell's equations, it is 
difficult to extend such techniques to include effects like 
the screening of intramolecular interactions by the elec- 
trodes in molecular junctions. 

In contrast to this, effective field theory (EFT) pro- 
vides a concise, systematic method of constructing a tt- 
electron Hamiltonian starting from first principles by per- 
forming an expansion in a small parameter and then im- 
posing symmetry constraints. The result contains a few 
physically meaningful parameters, which are then renor- 




FIG. 1. Two isosurfaces of the average 7r-electron density 
(tp'(x)ip(x)) depict the electronic structure of gas-phase ben- 
zene within 7T-EFT. 



malized to include the aggregate effect of the degrees of 
freedom not explicitly retained. In this article we pro- 
ceed along these lines, first expanding the full electronic 
Hamiltonian of a conjugated organic molecule in a ba- 
sis of atomic orbitals and then dropping terms involving 
energies far from the 7r-electron bond energy. 

Imposing symmetry constraints and performing an ex- 
pansion in powers of the interatomic bond length then 
allows us to construct an effective Hamiltonian for the ir- 
electrons in a conjugated organic molecule that accounts 
for the effects of cr-electrons virtually. 7 As an example 
of this, we consider the particular case of gas-phase ben- 
zene, for which we formulate an effective Hamiltonian 
with only four adjustable parameters: the on-site repul- 
sion U, the nearest-neighbor hopping matrix element t, a 
dielectric constant e, and the 7r-electron quadrupole mo- 
ment Q. In principle, these could then be renormalized 
ab initio, e.g. by using perturbation theory to freeze out 
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degrees of freedom far from the 7r-electron energy scale; 
however, since this is tedious and would not enhance the 
predictive power of our model, we instead fit the param- 
eters directly to experiment. 

Next, we show how screening from metallic electrodes 
can be incorporated into this scheme without introduc- 
ing additional parameters by considering the multipole 
moments of image charge distributions. We then use this 
method of screening to calculate the screened HOMO- 
LUMO gap of benzene near a metallic electrode, as well 
as formulate a realistic model of a gold-benzenedithiol- 
gold junction, including effects arising from the pres- 
ence of the thiol sidegroups. Finally, we also use this 
7r-electron effective field theory (71--EFT) to compute the 
7r-electron contribution to the van der Waals interaction 
between benzene and a metallic electrode. 



I. BARE HAMILTONIAN 

Using the Born-Oppenheimer approximation, the one- 
body term in the electronic Hamiltonian for an isolated 
molecule can be written as 

H{1) = E/ <Px^{x) (^v 2 + Vv(of) (l) 

where V is the interaction between the electrons and the 
atomic nuclei. The operator that creates an electron with 
spin a in the nth element of a basis of atomic orbitals 
{4> n } can be expressed as: 



4 CT = / d?x4> n mi{z) 



Multiplying this by the inverse of the overlap matrix 
S nm — (<Pn\<Pm) and summing over m implies 



m m 

= ^M (2) 

where we have made use of the completeness relation for 
a non-orthogonal basis 

4>n{^)S~^(j)* m (x) = 6(x-x l ) 

nm 

Combining equations ([T]) and © then gives: 

H ■ = ^ ^~(-nrndl l(7 d n i(T (3) 



where 



= J d 3 x <j>* n {x) V 2 + Vj <P m (x) (4) 



and 



(5) 



If we keep only nearest-neighbor terms this reduces to 
the tight-binding Hamiltonian 

H w = Ve d) d - V t d f d 

n {n,vn),<7 

where t nm = Hnm and e„ = Hnn . 

Similarly, the two-body term in the electronic Hamil- 
tonian can be written as 



H{2) = \Hl d^x^l{x x )^,{x 2 ) w 



X2\ 



which, in terms of the atomic orbital basis, is equivalent 
to 

H (2) = \ £ U nmkl dlJ m<J ,d ka ,d la (6) 



nraklaa' 



where 



U nmk l = J d Xid X 2 (/)* (Xl)0* n (f 2 ) 
<t>k{x2)4>l{xi) 



\Xl - X 2 \ 



and 



kl 



U nm kl — ^] SnoSmpU op qr{S kq ) (S [r ) . (7) 
opqr 

Together, equations ©-© and ©-© give the full 
electronic Hamiltonian from first principles: 

H = HM+H&\ 

but do so in terms of a basis that is impractically large 
for use within existing many-body techniques. To over- 
come this difficulty, in the next section we formulate an 
effective Hamiltonian in a reduced basis, explicitly re- 
taining only the degrees of freedom necessary to describe 
the 7r-electron dynamics. 



II. EFFECTIVE HAMILTONIAN 

The first step in constructing the effective Hamilto- 
nian is culling elements of the basis that lie far from the 
energy scale of interest. To this end, we first exclude 
atomic orbitals that do not participate in chemical bond- 
ing (those corresponding to core or excited electrons), 
which, in a 7r-electron system, leaves an effective s or- 
bital and three effective p orbitals at each atom. The 
former hybridize with the effective p x and p y orbitals 
giving rise to three sp 2 hybrids that form the a bonds 
between the atoms. The remaining effective p z orbitals, 
which, for a planar molecule, cannot hybridize with any 
of the cr-electrons without breaking inversion symmetry, 
are occupied by one electron on each atom and form 7r 



3 



bonds with weaker binding energies. Because of this en- 
ergy difference we also omit the atomic orbitals partici- 
pating in the a bonds, though this approximation could 
be relaxed at the expense of a larger basis. 

The effective Hamiltonian for the remaining effective 
p orbitals can then be determined using equations Q 
through ([7]) if the effective orbitals are known. In prin- 
ciple, these could be calculated directly, e.g. by using 
perturbation theory to freeze out the degrees of freedom 
far from the 7r-electron energy scale; however, as noted 
previously we find it more practical to parametrize these 
expressions by imposing symmetry constraints. 

To do this, we work initially in the asymptotic limit 
where the interatomic bond length is large compared 
to the size of the effective orbitals. This condition im- 
plies that matrix elements U n mki with n ^ I or m / I; 
and overlap integrals S nm with n ^ m are exponentially 
small, allowing us to reduce the interaction matrix (Eq. 
©) to 



Mnmkl — U nrn 



kl 



OnlOmk / d Xld X 2 



Pi - x 2 



(8) 



where <p n are now effective instead of bare orbitals. Al- 
though it is known^ that the terms neglected are not a 
priori negligible at typical interatomic distances, it has 
been suggested that this approximation can be justified 
by the use of orthogonalized orbitals. — Here we offer 
a simpler perspective more consistent with the spirit of 
EFT, namely that the neglected terms are accounted for 
virtually when the parameters in the Hamiltonian are 
renormalized. We also note that Eq. ([8]) is equivalent to 
the "neglect of differential overlap approximation" that 
has already been used extensively elsewhere, but that 
in the context of EFT it is simply the requirement that 
the effective Hamiltonian be local. It has been explicitly 
shown elsewhere- 2 that this is a satisfactory approxima- 
tion for 7r-conjugated systems. 

Expanding Eq. ([8]) in powers of the interatomic bond 
length yields a standard electrostatic multipole expan- 
sion, and, if we assume the effective p orbitals pos- 
sess azimuthal and inversion symmetry, U nm can be 
parametrized up to the quadrupole-quadrupole interac- 
tion in terms of the on-site repulsion U nn and the zz com- 
ponent of the quadrupole moment Q n associated with 
each orbital, as well as a dielectric constant e included to 
account for the polarizability of the a and core electrons. 

Explicitly, this gives 

Unm — U nn 8 nm 

+ (i — s nm ) KL M + u^ 1 + u%™ + ug£) 

+ 0(r- 6 ), (9) 

where JJ is the monopole-monopole interaction, JJ^ M 
is the quadrupole-monopole interaction, and U®Q is 



the quadrupole-quadrupole interaction. For two orbitals 
with arbitrary quadrupole moments and Qm sepa- 
rated by a displacement r, the expressions for these are 

„2 



u. 



MM 



cr 



jjQM \- Qij 



(10) 

(ii) 

(12) 
(13) 



ijkl 



where 

Wijki = SuSkj + SkiSij - 5r~ 2 (r k 5urj + r k nSij 

+ Skirjri + nSkjri + rkriSij) + 35r~ i r l r J rirk 

is a rank-four tensor that characterizes the interaction of 
two quadrupoles. 13 Altogether, this provides an expres- 
sion for the interaction energy that is correct up to fifth 
order in the interatomic distance. 

To further reduce the number of free parameters it is 
convenient to simplify the effective Hamiltonian by re- 
quiring it to satisfy particle-hole symmetry. Although 
this is not strictly necessary within the context of 7T-EFT, 
the success of Pariser-Parr-Pople type semi-empirical 
models-which implicitly assume particle-hole symmetry- 
suggests that it is a good approximation to do so. Taking 
this to be the case, Eq. (j4]) then gives for the one-body 
Hamiltonian: 

' -h 2 



H 



(i) 



d 3 x(j) n (x) 



2m 



V 2 +^(5?) ct> m (x), 



where Vi (x) is the effective potential due to the ionic hole 
at site I: 

J e\x-x'\ 
Using Eq. (J5J) then gives 

H {1) - 5 I e (at) - Vf/ ; I 4- CI - 8 )t 

V !*» / 

where we have defined the atomic on-site energy as: 

= J d 3 x4>* n (x) (^V 2 + Vnix^j <Mf). 

Defining p n = 5Z<t d na d n(T and rearranging the two-body 
term then yields: 

1 1 i — 7_j n I n /_ j L nm ( - l n ( 7 771(7 

n (n,m),(T 

+ U nm(Pn ~ l)(Pm - 1) 

nm 

~t~ ~ ^ ^ U nn p n — ^ ^ U nm . 
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Finally, adding the mutual repulsion of the ionic cores 
5 En/m fnm gives the full effective molecular Hamilto- 
nian: 

F-Ve (at) o - V t S d 

n (n,m),cr 



2 ^ ^ U nrn Q n Q 

nm 



where we have introduced the effective charge operator 
defined by q n = p n — 1. In conjunction with Eq. ([9]), 
this expresses the effective Hamiltonian for an arbitrary 
7r-electron system in terms of the tight-binding matrix 
t nm , the on-site repulsion U nn , a dielectric constant e, 
and the 7r-electron quadrupole moment Q n . 

In the remainder of this paper we focus on benzene 
as a benchmark system, in which case the Hamiltonian 
reduces to 



$ d 



2 ^ 

nm 



(14) 



where U nn = U and Q z n z = -Q* x /2 = -QHff/2 = Q by 
symmetry. The molecular chemical potential p is fixed 
by the experimental ionization energy^— and electron 
affinity ji^ 



IE- EA 



-4.06 eV, 



whereas the four other parameters must be renormalized 
by fitting to experiment, which is the subject of the fol- 
lowing section. 



III. RENORMALIZATION: FITTING THE 
GAS-PHASE SPECTRUM 

We have renormalized the parameters in our effec- 
tive Hamiltonian for gas-phase benzene by fitting to ex- 
perimental values that should be accurately reproduced 
within a 7r-electron only model. In particular, we have 
simultaneously optimized the theoretical predictions of 
1) the vertical ionization energy, 2) the vertical electron 
affinity, and 3) the six lowest singlet and triplet excita- 
tions of the neutral molecule. 

This was done by exactly diagonalizing Eq. (fT4")) with 
the interatomic bond length 2 ^ fixed at 1.40 A. In partic- 
ular, using the OQNLP algorithm 2 ^ for nonlinear global 
optimization we minimized the RMS relative error of our 
predictions for the quantities in the first column of Table 
|IJ The results of this procedure, which converged to the 
same solution regardless of initial conditions, are sum- 
marized in column two of the same table. The optimal 
parametrization for the 7T-EFT was found to be t — 2.70 
eV, U = 9.69 eV, Q = -0.65 eA 2 and e = 1.56 with a 
RMS relative error of 4.2 percent. 



Also appearing in TableQ]are the predictions of a recent 
Pariser-Parr-Pople type semi-empirical model£ as well as 
those of the original Ohno parametrization. 4,? Compared 
to the recent PPP model, 7T-EFT fits the optical spec- 
trum of gas-phase benzene to a similar degree of accu- 
racy and gives better results for the ionization energy 
and electron affinity. Moreover, the parameters common 
to both models have comparable values, namely those 
given above for our model and for the model of Castle- 
ton et al. t = 2.64 eV, U = 8.9 eV, and e = 1.28. 
The 7T-EFT on-site repulsion is also in qualitative agree- 
ment with recent RPA-based calculations of the effective 
Coulomb repulsion in graphene^ 

Although our effective quadrupole moment has no di- 
rect counterpart in phenomenological models, its value 
can be compared to the bare quadrupole moment of a 
hydrogenic 2p z orbital, which is given by 



-24e(a /Z) 2 



(15) 



where do is the Bohr radius and +Ze is the nuclear 
charge. Using this, we find that our 7r-orbital quadrupole 
moment corresponds to a hydrogenic p orbital bound by 
an effective charge of +3.22e, or, equivalently, with an 
effective Bohr radius of 0.16 A. This is consistent with 
the expectation that the sp 2 orbitals forming the a bonds 
provide only weak screening of the atomic core, which has 
a net charge of +4e. For the purpose of visualization, ef- 
fective hydrogenic orbitals can also be used to render the 
average 7r-electron density (tl)'(x)ijj(x)), as shown in Fig- 
ure [1] 



IV. SCREENING BY METALLIC 
ELECTRODES: THE IMAGE MULTIPOLE 
METHOD 

In this section, we extend the preceding model to in- 
clude the effect of screening by metallic electrodes, which 
for simplicity are modeled as planar or spherical conduc- 
tors. In the regime where the characteristic response time 
of the electrons in the electrode is much shorter than the 
timescale of the 7r-electron dynamics, this can be done us- 
ing the method of images via a straightforward extension 
of Eq. This is expected to be the case for conjugated 
organic molecules in the vicinity of gold electrodes, in 
which case the metallic plasma frequency 2 ^ oj p ~ 9 eV/H 
is large compared to the frequency scale of ir excitations 
0;^ R3 2t/h ?s 5 cV/h. The leading order correction to the 
metallic dielectric function, given by the GW approxima- 
tion, then goes as (w^/wp) 2 ~ 0.3. Explicit calculations 
using the GW approach also suggest that corrections to 
the image charge method tend to be small for organic 
molecules adsorbed on a metallic surfaced 

In the following subsections, the multipole moments 
of the image charge distribution generated by an orbital 
near planar and spherical conductors are described. To 
determine the screened interaction matrix, interactions 
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Exp.ii-.22 


tt-EFT 


PPP 

Castleton et. al^ 


PPP 
Ohnc ' 


Ionization Energy (eV) 


9.23 


9.26 


9.05 


9.78 


Electron Affinity (eV) 


-1.12 


-1.14 


-0.93 


-1.67 


1 N C U-LI eLl OJJcLLI Hill 1 c V 1 










Singlet 


4.90 


4.87 


4.76 


4.23 




6.21 


6.08 


6.30 


5.52 




6.93 


7.59 


6.93 


6.81 


Triplet 


3.93 


4.10 


3.99 


3.52 




4.75 


4.92 


4.74 


4.32 




5.60 


6.17 


5.84 


5.58 


RMS Relative Error (%) 


N/A 


4.2 


6.0 


19.0 



TABLE I. Experimental data for the vertical ionization energy,— — vertical electron affinity,— and optical spectrum,— — of 
gas-phase benzene compared to the predictions of 7T-EFT, a recent PPP-type model,— and the Ohno parametrization.— £ The 
best fit parameterization of our 7T-EFT was determined to be t — 2.70 eV, U = 9.69 eV, Q = —0.65 eA 2 and e = 1.56. 



between these and the orbital multipole moments are in- 
cluded in U nm using equations (ITU1) through (flU)) . Over- 
all, the two-body Hamiltonian should give the energy re- 
quired to prepare the molecular charge distribution by 
bringing each of the electrons in from infinity with the 
electrodes maintained at fixed electrostatic potentials. 
This can be ensured using a number of different count- 
ing schemes, but we take one that ensures a symmetric 
interaction matrix, namely 

Unm = U nm + 6 nm U n ^ + — (1 — 5 nm )(U n ^ l + U mn ), 

where U nm is the unscreened interaction matrix, Unm is 
the interaction between the nth orbital and the image 
of the mth orbital, and U nm is the screened interaction 
matrix. Since the image multipole moments of an or- 
bital change as it is brought in from infinity, one might 
expect a pref actor of 1/2 in the second term of the pre- 
ceding equation, however, this is already present in the 
Hamiltonian itself. 

When multiple electrodes are present, the image of an 
orbital in one conductor produces images in the other 
electrodes, resulting in an effect reminiscent of a hall of 
mirrors. We deal with this by including these "higher 
order" multipole moments iteratively until the difference 
between successive approximations of U nm drops below 
a predetermined threshold. In practice, this procedure 
converged rapidly. 

Within the foregoing scheme, the case where one or 
more electrodes are maintained at a fixed potential other 
than zero can be treated straightforwardly by including 
image charges that contribute to the one-body Hamilto- 
nian rather than to U nm . For example, a spherical con- 
tact with radius R at potential V can be treated using 
a hypothetical point charge q = VR at the center of the 
electrode. This technique is especially useful for trans- 
port calculations in the context of molecular junctions, 
as it provides the full junction Hamiltonian at finite bias, 



alleviating the need for the phenomenological models of 
capacitive lead-molecule coupling that have been relied 
upon on in the pastil 



A. Screening by a planar electrode 

In classical electrostatics, the image of a charge dis- 
tribution near a planar conductor is merely the mir- 
ror image of the charge distribution itself. Thus an 
orbital with monopole moment q and quadrupole mo- 
ment Q lJ located a distance r away from a conducting 
plane produces an image orbital inside the conductor lo- 
cated at depth r with multipole moments q — —q and 
Qij — —^ kl TikTjiQ kl , where Tifc is a transformation 
matrix representing a reflection about a plane parallel to 
the surface of the conductor, i.e. 

Tik = Sik - 2hihk, (16) 
where n is the unit vector normal to the planar surface. 



B. Screening by a spherical electrode 

An orbital with monopole moment q and quadrupole 
moment Q 13 located a distance r from the center of a 
spherical electrode with radius R induces an image dis- 
tribution at f = R 2 /r with monopole and quadrupole 
moments 

r Ir L — ' 

ij 

and 




respectively, where is a transformation matrix repre- 
senting a reflection about the plane normal to the vector 
f, similar to Eq. (|f 6p . 

Thus the orbital quadrupole moment induces a higher 
order image monopole moment, as well as an image of 
itself that is deformed and reflected. An image dipole is 
also generated, but its interaction with the orbital charge 



distribution is of order 
here. 



and so we have neglected it 



V. SCREENING OF THE HOMO-LUMO GAP 

Although 7T-EFT could be used to study a wide variety 
of phenomena involving conjugated organic molecules, 
our primary motivation in formulating it has been to fa- 
cilitate realistic many-body calculations of transport phe- 
nomena in molecular junctions. In particular, while re- 
cent semi-empirical models^ reproduce the low-lying ex- 
citations of gas-phase benzene, their predictions of quan- 
tities relevant to transport, namely the HOMO-LUMO 
gap and the optical excitations of the ionized molecule, 
are less accurate. Moreover, in a molecular junction these 
quantities are renormalized by screening from metallic 
electrodes as well as the presence of linker groups not 
explicitly included in the molecular Hilbert space. 

Within 7T-EFT these effects can be clearly seen: Con- 
sider the spectral function of gas-phase benzene, which 
we evaluate at the many-body level using the non- 
equilibrium Green's function formalism as described in 
appendix A. Figure [5] shows this quantity, along with ex- 
perimental values for the vertical ionization energy (9.23 
eV), vertical electron affinity (—1.12 eV), and the first 
optical excitation of the cation (3.04 eV). As a guide to 
the eye, the spectrum has been broadened artificially us- 
ing a broadening matrix of T nm = (0.2 eV)S nm . As an 
aside, we note here that the close agreement between the 
experimental values and the maxima of the spectral func- 
tion suggests our model is accurate at this energy scale. 
In particular, the accuracy of the theoretical value for 
the lowest optical excitation of the cation is noteworthy, 
as this quantity was not fit during the renormalization 
procedure but rather represents a prediction of 7T-EFT. 

Screening effects become evident when the molecule is 
brought into proximity with the surface of a planar elec- 
trode. Figure [3] shows the reduction of the ionization 
energy and electron affinity as a function of electrode- 
molecule distance in this scenario, and the HOMO- 
LUMO gap, given by IE — EA, is reduced commen- 
surately. These results, based on the image multipole 
method, are also consistent with recent GW-based inves- 
tigations of screening! 28 ' 30 

We also considered the prototypical benzene-gold junc- 
tion, consisting of benzene linked to two gold electrodes 
via thiol side groups. Although this junction can oc- 
cur with a wide variety of different geometries, in this 
example we have taken the configuration shown in Fig- 
ure [31 The electrodes are modeled as metallic spheres 
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FIG. 2. The spectral function of gas-phase benzene broadened 
artificially as a guide to the eye. The dashed orange lines are 
fixed by (left to right) the lowest lying optical excitation of 
the molecular cation,— -18 ' 31 the vertical ionization energy of 
the neutral molecule,—^ and the vertical electron affinity of 
the neutral molecule. 19 
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FIG. 3. The ionization energy and electron affinity of benzene 
oriented parallel to the surface of a screening plane, shown as 
a function of distance. 



with radii of 0.5 nm, and the partially ionic character 
of the gold-sulfur bond has been accounted for by plac- 
ing point charges of — 0.67e at the locations of the sulfur 
atoms. The latter value was determined in conjunction 
with the tunneling-width matrix (rn = = 0.44 eV) 
via a simultaneous fit of the experimental thcrmopowcr 32 
and conductance, 33 using the techniques described in ap- 
pendix A. The upper panel of Figure [5] shows the spec- 
tral function for this junction in the simple case where 
the tunneling-width matrix is the same as in Figure [5J a 
choice which simplifies comparison of the two cases. 
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FIG. 4. The geometr y 3 ' for the benzenedithiol junction as- 
sociated with the spectral function shown in Figure [5] The 
electrodes have been placed so that the screening surface 
lies one covalent radius^ beyond the position of the outer- 
most gold nuclei, a convention that has been investigated 
elsewhere 3 - in the context of atom-surface van der Waals in- 
teractions. 



Screening from the electrodes reduces the HOMO- 
LUMO gap by 12.5 percent as compared to the gas-phase, 
and the dipole formed by the gold-sulfur bond shifts the 
chemical potential of the molecule up by 1.4 eV. For com- 
parison, we have also calculated the spectral function 
of the same junction, but with the electrodes modeled 
as planes (Figure [SJ bottom) , in which case the screen- 
ing is maximal and the HOMO-LUMO gap is reduced 
by 19 percent. These results are qualitatively consistent 
with recent GW-based investigations of screening effects 
wherein a molecule is adsorbed on a metallic surface ) 28 i 30 
though the scenario we have considered is not directly 
comparable to these works. For both electrode geome- 
tries a splitting of the two-fold degenerate HOMO and 
LUMO resonances can also be seen, which arises from 
the interaction between the 7r-electrons and the dipoles 
associated with the partly ionic gold-sulfur bonds. 

VI. tt-ELECTRON CONTRIBUTION TO THE 
VAN DER WAALS INTERACTION 
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Junction (R = 0.5 nm) 
Gas phase 




Junction (R = °o) 
Gas phase 




E(eV) 

FIG. 5. Top: The spectral function of the Au-1,4- 
benzenedithiol-Au junction depicted in Figure[3]at room tem- 
perature, together with the gas-phase density of states from 
Figure [21 To facilitate comparison, the same broadening has 
been used in both cases. The dotted orange line at —5.1 eV 
indicates the position of the experimental chemical potential 
of clean gold.— Bottom: The spectral function of the same 
junction with planar instead of spherical electrodes. 



As a second application of the image-multipole 
method, we consider the 7r-electron contribution to the 
van der Waals interaction between a molecule and a 
metallic electrode. Experimentally, such interactions are 
important when a molecule is adsorbed on a metal sur- 
face, or in single-molecule junctions in which a molecule 
bonds directly to metallic electrodes, as in the Pt- 
benzene-Pt junctions investigated recently by Kiguchi et 
a \39. Theoretically, the van der Waals interaction also 
represents a unique challenge in that it is a true many- 
body phenomenon arising from quantum correlations in- 
duced by long-range interactions. As such, it is outside 
the scope of local approximations to density functional 
theory, and modeling van der Waals interactions using 
nonlocal functionals is a topic of ongoing research.— ~— 
In contrast to this, the preceding treatment of screening, 
in conjunction with a full many-body treatment of the tt- 
electrons on the molecule, makes it possible to calculate 
the 7r-electron contribution to van der Waals interaction 
straightforwardly with no extra adjustable parameters. 

In particular, by exactly diagonalizing the few-body 
molecular Hamiltonian with and without the effects of 
screening included in U nm , it is possible to infer the van 
der Waals interaction at arbitrary temperature between 
a molecule and a metallic electrode by comparing the 
expectation values of the Hamiltonian in these two cases: 

E v dw = (H) — (H) 

This procedure was carried out at zero temperature 
for benzene oriented parallel to the surface of a spherical 
electrode over a large range of electrode-molecule dis- 
tances, and the results are shown in Figure [5] When 
the molecule is near the surface of the electrode E v dw = 
— which is the expected asymptotic dependence for 
the van der Waals interaction between a molecule and 
a planar conductor. Conversely, when the molecule is 
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FIG. 6. The 7r-electron contribution to the van der Waals 
interaction between benzene and a spherical electrode with a 
radius of 10 nm, plotted as a function of the distance from 
the conducting surface. At all distances the molecule is ori- 
ented parallel to the surface of the electrode. The dashed 
green and orange lines show the expected asymptotic depen- 
dence in the near and far fields respectively. Inset: The same 
quantity very near the surface of the electrode, including a 
phenomenological 43 gold-carbon hard-core repulsion. 
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FIG. 7. The orientation dependence of the 7r-electron contri- 
bution to van der Waals interaction between a planar elec- 
trode and a benzene molecule centered 2 nm from the metal 
surface. The molecule is initially oriented parallel to the elec- 
trode and then rotated by an angle 9±_ about the axis perpen- 
dicular to the plane of the molecule, followed by a rotation of 
fly about an axis within the plane of the molecule. 



C3 ps 1.56 eV A 3 . We also investigated the orientation 
dependence of the van der Waals interaction between a 
planar electrode and a benzene molecule, as depicted in 
Figure [71 which shows a significantly stronger attractive 
interaction when the plane of the molecule is oriented 
perpendicular to the surface of the electrode. 

The van der Waals coefficient C3 is fundamentally re- 
lated to the molecular polarizability tensor onj. Thus, for 
an axially symmetric molecule such as benzene, a simpli- 
fied single-oscillator model can be used to derive semi- 
empirical formulae relating ay to C3 with the molecule 
oriented either parallel or perpendicular to the surface of 
a planar electrode^ 



ct 



Ed 
32 
Ed 
32 



(2au + 2ay) 
(a± + 3cif|| ) 



(17) 
(18) 



Here Ed is the energy of the principal dipole-allowed 
optical transition, and an and a± are respectively the 
molecular polarizabilities parallel and perpendicular to 
its plane of symmetry. Calculated within 7T-EFT, Ed = 
7.59 eV, ay = 3.24 A 3 and a± = 0.00 A 3 . These fig- 
ures in conjunction with the preceding expressions imply 
Cj ps 1.54eVl 3 and C^/c\ = 1.5, which agree very 
closely with the values obtained by direct calculation us- 
ing tt-EFT. 

Experimentally, an = 12.31 A 3 , a± = 6.35 A 3 and 
Ed = 6.93 eV for benzene, and in this case Eq. (JT7J) gives 
Cj «8.08 eV A 3 , which is roughly five times larger than 
that predicted by 7T-EFT. This discrepancy can be at- 
tributed to the significant contribution of the cr-electrons 
to the molecular polarizability, as evidenced by the large 
experimental value of a±, which arises from a-ir transi- 
tions. This underscores the importance of the tr-electron 
dynamics in the context of van der Waals interactions, 
which arise from the long-range spatial correlation of 
purely virtual processes. In contrast to this, the effect 
of the a-electrons on real 7T-7T transitions, such as those 
involved in transport, should be well described by 7T-EFT. 
Moreover, as noted previously, the cr-electrons dynamics 
can be explicitly included within effective field theory at 
the expense of a larger Hilbert space, and we believe that 
such a 7TC7-EFT would accurately reproduce the full van 
der Waals interaction between a conjugated molecule and 
a metallic electrode. 



VII. CONCLUSIONS 



far from the electrode E v dw = — which is the usual 
asymptotic dependence of the van der Waals interaction 
as given by the Lennard- Jones potential. A clear tran- 
sition between the two regimes can be seen around 10 
nm, the radius of the electrode. In the near-field region 
the constant of proportionality predicted by 7T-EFT is 



We have shown how EFT can be used to provide a con- 
cise derivation of an effective Hamiltonian for 7r-electron 
systems by performing a multipole expansion, imposing 
symmetry constraints, and then renormalizing a few ad- 
justable parameters. In particular, we have optimized 
the parameters appearing in an effective Hamiltonian for 
gas-phase benzene, Eq. (|14p . by fitting to experimental 
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data for 1) the vertical ionization energy, 2) the vertical 
electron affinity, and 3) the six lowest singlet and triplet 
excitations of the neutral molecule. This procedure yields 
a fit which is comparable to or better than traditional 
PPP models^i£, and gives U — 9.69 eV for the on-site 
repulsion, t — 2.70 eV for the nearest-neighbor hopping 
matrix element, e = 1.56 eV for the dielectric constant, 
and Q = —0.65 eA 2 for the 7r-electron quadrupole mo- 
ment. These values of U, t, and e are consistent with 
those used in previous 7r-electron models£&2£, while Q is 
a new physical parameter in our approach, which takes 
the place of the ad-hoc functional forms assumed in PPP 
models and governs the corrections to 1/r interactions at 
short distances. 

We have also utilized 7T-EFT to model the screening 
of intramolecular Coulomb interactions by nearby metal- 
lic electrodes. Within our approach, lead-molecule cou- 
pling is treated using a two-step process wherein all long 
range Coulomb interactions are included nonperturba- 
tively before lead-molecule tunneling is accounted for via 
Dyson's equation. The ability to include screening ef- 
fects via image multipoles-without additional adjustable 
parameters-represents a significant advantage of 7T-EFT 
over PPP models, which utilize interactions that do not 
satisfy Maxwell's equations. 

In particular, we have shown how 7T-EFT facili- 
tates a realistic description of the prototypical Au-1,4- 
benzencdithiol-Au junction. In this case, we obtain es- 
timates of both the HOMO-LUMO gap and the chem- 
ical potential of the molecule in the junction. Beyond 
this, we have also described how 7T-EFT facilitates trans- 
port calculations at finite bias without the introduction 
of additional adjustable parameters, which we foresee as 
a fruitful area for future work. 
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Appendix A: Many-body theory of transport in 
molecular junctions 

Within the non-equilibrium Green's function approach 
to studying transport in molecular junctions, a quantity 
of central importance is the retarded Green's function G 
of the molecule coupled to the electrodes. In the energy 
domain and using matrix notation, this can be expressed 
via the Dyson equation as: 



G — G mo i + G mo /SG, 



(Al) 



where G mo i is the interacting Green's function of the 
molecule without tunnel coupling to the electrodes, but 
including long-range Coulomb interactions between the 
7r-electrons and their image multipole moments in the 



leads. The self-energy £ can be partitioned into the tun- 
neling self-energy Et associated with the lead-molecule 
bonds, and a correction to the Coulomb self-energy AEc 
arising from lead-molecule coherence: 



Far from resonance and at room temperature AEc ~ 
0, and so in the present context we neglect this 
correction-an approximation which is justified in detail 
in ref. [2^. Assuming the leads can be modeled as Fermi 
liquids with good screening, the electron-electron inter- 
actions within them can be neglected and the tunneling 
self-energy associated with a given electrode can be ex- 
pressed as 

E T = Vg{E)V\ 

where g(E) is the retarded Green's function of the lead 
and V n k are the matrix elements coupling the lead and 
molecule. In the broad-band limit wherein the density 
of states in the electrodes varies slowly in the vicinity of 
the metallic Fermi level, the self-energy then reduces to 
a purely imaginary matrix with no energy dependence^ 



(A2) 



Here the tunneling-width matrix r Q associated with lead 
a given is by 



r 



na.mcr 



2n P {e f )V n V^5 a 



where p{ei) is the density of states at the metallic Fermi 
level, and V n is the matrix element between the nth n- 
orbital in the molecule and the lead states in the vicinity 
of the Fermi level. The diagonal elements of this equa- 
tion are equivalent to Fermi's golden rule, with Tr {T a /h} 
giving the rate at which electrons in lead a are being in- 
jected into the molecule. 

Aside from the self-energy, the other ingredient needed 
to evaluate Eq. (|Alj) is the Green's function of the iso- 
lated molecule. This is determined exactly by first find- 
ing the few-body eigenstates and eigenenergies E v 
of the gas-phase molecule, and then using these to ex- 
plicitly evaluate the molecular Green's function : 29 ! 46 



G n 



E 



[P(u)+P(v')]C(v, v') 
E - (E u > - E v ) + M+ 



(A3) 



Here P(v) is the statistical occupancy of the uth. eigen- 
state, given at equilibrium by the grand canonical ensem- 
ble, and 

C na ,m<j'{v,v') = {v\d n <j\v'){v'\d ] m(jl \v) 

are many-body matrix elements, where, in the present 
context, d\ na creates an electron with spin a in the mth 
7r-orbital of the molecule. 
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Altogether, equations (|A1[) . (|A2[) and (|A3[) provide a 
method for obtaining the full interacting Green's function 
of the molecule coupled to the electrodes, which may then 
be used to calculate the various physical quantities of 
interest. For example, the spectral function is given by>^> 



A(E) = -2 ImG, 



ular Green's function via the expression: 

T a f3 — Tr {T a GT pG^} , 

where T a and are the tunncling-width matrices as- 
sociated with leads a and /3 respectively. This quantity 
may then be used to evaluate the various electronic trans- 
port quantities of interest,— such as the elastic electrical 
current 



the trace of which is proportional to the effective single- 
particle density of states: 



p(E) = —Tr{A} 



and elastic thermal current 

# = \ E / dE ( E - f i °) T °p (fp - f") 



Similarly, the elastic transmission function between 
two electrodes can also be obtained from the full molec- 



flowing into lead a. Here f a (E) and fi a are respectively 
the Fermi-Dirac distribution and chemical potential as- 
sociated with lead a. 
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